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Ph Abstract 

< 

■^ We consider operators in the domains with the boundaries and 

J« derive sharp spectral asymptotics (containing non-Weyl correction) 

in the case when Hamiltonian flow is periodic. 

Even if operator is scalar but not second order (or even second- 
order but there is an inner boundary with both refraction and reflec- 
tion present) Hamiltonian flow is branching on the boundary and the 
^T) notion of periodicity becomes more complicated. 
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Q Introduction 



IL" In section 6.2 of [Ivr2]^'' 

•^ we considered the case when there is no boundary and Hamiltonian 

rS flow is periodic and derive very sharp spectral asymptotics. Now we want 

ci to achieve similar results when there is a boundary. We do not expect 

asymptotics to be that sharp, but better than 0{h^~'^) and have non-Weyl 

correction term. However the presence of the boundary, more precisely, 

branching of the rays on the boundary brings the new possibilities. 



^) This article is a rather small part of the huge project to write a book and is just 
section 8.3 consisting entirely of newly researched results, section 6.2 roughly corresponds 
to Section 4.7 of its predecessor V. Ivrii [Ivrl]. External references by default are to 

[Ivr2]. 
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Even if operator is scalar but not second order (or even second-order 
but there is an inner boundary with both refraction and reflection present) 
Hamiltonian flow is branching on the boundary and the notion of periodicity 
becomes more complicated. 

For simplicity we consider only Schrodinger operators and derive sharp 
spectral asymptotics (containing non-Weyl correction) in the case when 
Hamiltonian flow is periodic. 



1 Discussion and plan 

We start our analysis in subsection 2 from the case when there is no branching 
and the Hamiltonian flow (with reflections) is periodic. Then if period is Tq, 
M^j-g = / we have the same equality 

(1.1) e'^'''^°^Q = e''^°^Q 

as before where supp q is disjoint from the boundary and only transversal to 
the boundary trajectories originate at suppg and L is /?-pseudo-differential 
operator. Then 

(1.2) F,^,-i.x(^ - n)r{Q,,uV2y) = {2T^hf-'j{n) + 0{h^-'+') 



T, 







with 



(1.3) ^^"^ = y e'"'q,q2d^r 

as X is supported in (— |, |) and equals 1 in ( — |, |) and n E Z, \n\ < h^^. 
Here as before E,- = {a(x, ^) = r} and (x^ = dxd^ : da is the measure there, 
£ is the principal symbol of L. 

If i is "variable enough" on Z^ so the right-hand expression is decaying 
as n — 7- oo we can improve remainder estimate 0(/)"^^'^) (to what degree 
depends on the rate of the decay). In section 6.2 essentially i was defined by 
integral of the subprincipal symbol along periodic trajectories. Now however 
i can pick up extra terms at the moment of reflection: for example, for 
Schrodinger operator increment of £ is and vr for Neumann and Dirichlet 
boundary condition respectively; however more general boundary condition 
brings increment which depends on the point {x',^') G 7"*/ of reflection. 

Then in section 3 we consider a more complicated case when two operators 
are intertwined through boundary condition and all Hamiltonian trajectories 
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of one of them are periodic but of another is not and, moreover, if the generic 
Hamiltonian bilhard is periodic, it is in fact Hamiltonian billiard of the first 
operator. Then (1.1) fails but (1.2)-(1.3) is replaced by (1.2), (1.4) 

(1.4) J{n) = [ e'"'-\"\''q,q2dyi^,^ 

where Zi ^ and yiir correspond to the first operator and £' > depends 
on the portion of energy which was whisked away along trajectory of the 
second operator at reflection point. Then if 

(1.5) m,r({(x,0 G Zl,r,^' = 0, V£ = 0}) = 

we conclude that J{n) = o(l) as n — )> oo and we will be able to improve 
remainder estimate 0(/7-'-~'') to o[h^^'^) and again there will be a non-Weyl 
correction term but it will depend also on i and i'. 

Finally, in subsection 3 we consider even more complicated when two 
operators are intertwined through boundary condition and Hamiltonian 
trajectories of both of them are periodic and moreover after a while all 
trajectories issued from some point at T* Y assemble there. Then we can 
construct a matrix symbol and its eigenvalues play a role of symbol i. 

2 Simple Hamiltonian flow 
2.1 Inner asymptotics 

Let U be an open connected subset in T*X disjoint from dX. Consider 
Hamiltonian billiards issued from U, let •t'f denote generalized Hamiltonian 

flow. Assume that 

(2.1) |Vx,ca(x,OI > e Vt V(x,0 G ct),(W), 

(2.2) |0(x)| + |{a,0}(x,OI > e Vt V(x,0 G MU), 

(2.3) cD,(x,0 = (x,0 with t = t(x, > V(x, G W, 
where G "^"^(X), > in X and (f) = on dX and 

(2.4) Along <t>t(W) reflections are simple without branching, i.e. L~^L{y,ri) fl 
{a{y, 7]) = t} consists of two (disjoint) points (y, rj^) such that ±{a, 0} > 
as (y, 7]) e ^t(W), y G dX and r = a(y, rj). 
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Let us recall that l : T*X\gx -^ T*dX is a natural map. 
Then exactly as in the case without boundary 

(2.5) One can select t{x, ^) = T{a{x, ^)) in U with T G ^°° (so period 
depends on energy level only^^ 

and 

(2.6) As (x, e W and < t < T(x, ^ ^t{x, meats 9X exactly (A/ - 1) 
times with N = const; so as A/ > 1 billiard consists of N segments with the 
ends on the boundary. 

We can consider instead of A with domain 

(2.7) ^{A) = {ueJ€"^{X),dBu = 0} 
operator f{A) with 

(2.8) f{T) = JT-\T)dT 

and introduce propagator e'^ ^^^^"^ so that the Hamiltonian flow M^f of f{A) 
is periodic with period 1. 

However now structure of f{A^ near dX is a bit murky; but it is not 
really serious obstacle since we can consider problem 

(2.9) {f{hDt)-A)e'^'''^^^^=Q, 

(2.10) gee'''"''^(^) = 

where f is an inverse function. Therefore its Schwartz kernel U{x,y,t) 
satisfies problem (2.9)-(2.10) with respect to x and transposed problem with 
respect to y. 

Note that due to assumptions (2.1)-(2.2) e'^ ^^^^^Q is a Fourier integral 
operator with the symplectomorphism v|/f provided Q is /)-pseudo-differential 
operator with the symbol supported in U. Further, due to (2.3), (2.8) M/i = / 
and therefore e'^ ^^^^Q is /?-pseudo-differential operator. So, we arrive to 



^^ But there could be subperiodic trajectories with periods n ^T(a(x,^)), n — 2,2,, 
The structure of subperiodic trajectories described in subsection 6.2.5. 
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Proposition 2.1. Let A be a Schrodinger operator with the scalar principal 
symbol a[x , ^) satisfying (2.1)-(2.3) and (2.8). Then 

(2.11) ^ih-H{AyiK ^ ^ih-hL ^^ ^ 

with e = h and h-pseudo-differential operator L; here n is Maslov's constant 
and the principal symbol i of L is defined by 



(2.12) e'-^= n 



e e '' 



l<k<N 

where i^ corresponds to k-th segment (see theorem 6.2.3) and i[ corresponds 
to k-th reflection. 

Example 2.2. Consider Schrodinger operator. Assume that at the reflection 
point a(x, = e? + a(x', C'), ^ = 6 + iPiO- Then 

(2.13) e'''^ = -{{r- a')^ - iP)~\{r - a')5 + i/3) 

with the right-hand expression calculated in the corresponding reflection 
point. 

Taking /3 = or /3 = oo (formally) we get i'l^ = and i'l^ = vr for Neumann 
and Dirichlet boundary conditions respectively. 

Then we can use all the local results of section 6.2, assuming that (2.11) 
holds with h' < e < h and derive local spectral asymptotics with the same 
precision (up to 0{h^^'^)) as there. 

Corollary 2.3. Let conditions of proposition 2.1 be fulfilled. Then asymp- 
totics with the remainder o[h^^^) holds provided 

(2.14) Mii^'O e r., Vz^£(x,o = 0}) = 

where Vj:^ means a gradient along Z^; 

(ii) Asymptotics with the remainder 0{h^~^^^) with small exponent S > 
holds provided 

(2.15) [ir{{{x.O e ^r, |V7(x,0l < e}) = oie'') 

with the small exponent 6' > 0; 
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(in) Furthermore asymptotics with the remainder o[h^~'^) holds as (2.1)-(2.2) 
are replaced by their non-uniform versions 

(2.2)' |0(x)| + |{a,0}(x,OI>O Vt W{x,^)eMUnZr\K), 

where |J.T(Ar) = 0. 

We leave to the reader to formulate statement similar to (iii) but with 
the remainder estimate 0{h^^^^^). 

Usually our purpose is the asymptotics with Qix = 'iplix), Qiy = "ipiiv) 
where V'y are smooth functions but not vanishing near the dX. 

Corollary 2.4. Asymptotics with the remainder estimate o{h^~'^) holds 
provided conditions (2.1), (2.2) , (2.3) and (2.14) are fulfilled inU = T*X, 
Qj = i/jj G ^°^(X) compactly supported. 

Proof. Due to corollary 2.3(iii) contribution of the zone {xi > e} to the 
remainder is Oi;[h^^'^) while contribution of zone {xi < e} to the remainder 
does not exceed 0{eh^^'^) as e > h. D 

These arguments could be improved to e = h^ and the remainder estimate 
could be improved to 0{h^''^^^). We leave the precise statements and 
arguments to the reader. Our goal is to derive more sharp remainder 
estimate. 

2.2 Asymptotics near boundary 

Let us consider X = {xi > 0} and U open subset in 7*]R^r"^ x [0, £:]xi x M^ 
where e > is a very small constant. We are interested in 

(2.16) F,^,-^,XT{t)Jr{e''-'''^^^Q) 

with Q = Q{x, hD'^, hDt) with symbol supported in U; let us rewrite it as 

with Q' = e'^'''f(^)Qe-''^~''f^^\ 
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Note that if (2.1), (2.2) are fulfilled, e > and t > Ce are small enough 
then Q' = Q'{x, hD^, hDt) is /)-pseudo-differential operator with the symbol 
supported in {e't < Xi < C't} and we can rewrite 

Q' = Q"{x, hD^) + Q"\x, hD^, hDt){hDt - f{A)); 

therefore we can rewrite (2.16) in the same form but with Q replaced by Q". 

Now we can apply all the arguments of the previous subsection and of 
6.2 and derive asymptotics with the remainder estimate as sharp as 0{h^^^) 
provided conditions (2.1)-(2.3) are fulfilled in the full measure as well as 
corresponding conditions of section 6.2 to i. Here as usual we lift the points 
of (/, r]', t) e T*dX to (/, T]^) e E,. 

So, we need to analyze only near tangent zone {xi < e, \{a, xi}\ < e}. 

Unfortunately the theory of the manifolds with the boundaries and 
periodic Hamiltonian (or even geodesic) flow with reflections is completely 
undeveloped. 

Let us assume that (2.1), (2.2)' and (2.3) are fulfilled for all (x, ^ e T*X : 
Ti < a(x,^) < T2 with Ti < T2. Then trajectories tangent to dX cannot 
penetrate into X and therefore xi = 0, ri < a'(x, ^') < T2 is incompatible 
with a;jx,0 <0. 

Further, xi = 0, ri < a'(x, ^') < T2 is incompatible with a'^^{x, ,^) > as 
well; really otherwise almost tangent to boundary billiards make very small 
jumps which contradicts to periodicity after N reflections with flxed N. 

Therefore dX is bicharacteristically flat 

(2.17) xi = 0,n< a'{x,C) <T2 ^ a;^(x,0 = 0. 
Further, (2.3) implies that 

(2.18) AH solutions of equation 

(2.19) z" + 6(vl/,(x',Ok = with 6(x',O = ^<x,(0,x',r) 
are To-periodic for any (x',^') G T*dX, ti < a'(x,^') < T2. 
Therefore (2.3) implies that 

(2.20) p(x, ^) o M/f X p with p = xi + |{a, xi}| as n < a{x,^) < T2. 
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Consider zone Ue = {{x,C) '■ p(-*^. ^ ^}- Blowing up (xi,,^i) — )■ 
£^^{xi, ^i) we can prove easily that after this 

(2.21) vl/,(x, = (vl/;,,(x, 0, ^l,^',t)Axi' 6)) , 
with 

(2.22) vl/;,,(x, ~ Koix'. a + E ^"Kin){x. 0. 

(2.23) %,e,)Axi'^i) ~ E^'''^(x',f,t),w(^i'^i) 

where ^[q is the Hamiltonian flow on T*dX and ^t n) is the linearized with 
respect to (xi,,^i) billiard flow (described by (2.19)). One can prove easily 
that U'j ,. and U''/ , ^, ^. . ^ are uniformly smooth. 

Then in blown- up coordinates (2.11) holds provided e > h2^^ because 
h = e~^h and then we can apply all the previous arguments and derive 
asymptotics; the remainder estimate is as sharp as 0{h^~^) provided (in not 
blown-up coordinates) 

(2.24) |V(x,c)^(x, 01 > eo as xi x p(x, < ei; 

X| X p{x,^) means exactly that xi > e|.^i|. 
Therefore 

Proposition 2.5. Let A be a Schrodinger operator with the principal sym- 
bol 3{x,^) satisfying (2.1) -(2.3) and (2. 17) -(2. 18). Let the subprincipal 
symbol and the boundary condition be such that (2.24) is fulfilled. All these 
conditions are supposed to be fulfilled in {(x,,^) : p{x,C,) < eo}- 

Then the contribution of zone {hi^^ < p{x,^) < 62} to the remainder is 

Since the contribution of zone {p(x, C,) < s} to the remainder does not 
exceed Ce^h^~^ and we can take e = h2~^ we conclude that the contribution 
of zone {p{x,^) < €2} to the remainder is 0(/)^~'^~^). One can also estimate 
this way the contribution of subperiodic trajectories. Thus we arrive (details 
are left to the reader) to 

Theorem 2.6. Let A be a Schrodinger operator with the principal symbol 
3{x,C,) satisfying (2.1) -(2.3) and (2. 17) -(2. 18). Let the subprincipal sym- 
bol and the boundary condition be such that (2.24) is fulfilled. All these 
conditions are supposed to be fulfilled in {(x,,^) : p(x, ,^) < eo}. 
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Then contribution of zone {p(x, ^) < £2} to the remainder is 0(/)^^^^^) 
while its contribution to the principal part of the asymptotics is given by 
the standard two term Weyl formula plus a correction term constructed 
according to section 5.2 for symbol L 

Problem 2.1. Recover 0{h^~'^) remainder estimate. This should not be ex- 
tremely hard especially if there are no subperiodic trajectories but definitely 
worth of publishing. 

Example 2.8. Consider Laplacian /?^A on the standard hemisphere §^ := 
{|x| = l,xi > 0} in W^^'^ or harmonic oscillator /)^A + |xp on the standard 
half-space M^ = {xi > 0}. 

(i) Consider boundary operator of example 2.2. Then A/ = 2 and £(z, r) = 
£'(z, r) -|- £'(f)(z), r) where z = (x', ^'), q[z) is the antipodal point and d' is 
defined by (2.13). Then one can express conditions to (. as conditions to /3; 
we leave exact statements to the reader; 

(ii) Consider Dirichlet or Neumann boundary conditions. Then as pertur- 
bation is hxi one can calculate easily that i = up + 0{p^) where p is the 
incidence angle of trajectory and n > 0. Then condition (2.24) is fulfilled. 

2.3 Discussion 

As we mentioned almost nothing is known about manifolds with all billiards 
closed; we are discussing geodesic billiards but one can ask the same questions 
about Hamiltonian billiards or Hamiltonian billiards associated with the 
Schrodinger operator; in the two last cases phase space S*X is replaced by 
the portion of the phase space {ti < a(x, ^) < r2}. We assume that the 
phase space is connected: 

Problem 2.9. (i) So far our only examples are a standard hemisphere §^ 
and a standard half-space M^ (see figure below); are any other there really 
different examples? 

(ii) In particular, are there manifolds with A/ 7^ 2 (i.e. with N = 1,3, 4, ... )? 
We can consider N copies Xi, ... , Xa/ of the standard quarter- sphere §^ fl {0 < 
^ < f } n {0 < < tt}, then glue {X/,,0 = vr} with {X^+i.^ = 0} for 
k = 1, ... , N, where X/v+i := X/y. Then for resulting manifold we have N 
refiections but for A/ 7^ 2 there will be singularity at North Pole. 
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(a) Hemisphere, projected to equatorial plane (b) half-plane 

Figure 1: Different billiards (actually, their x- projections) are sliown by 
different colors, tliey reflect from the boundary at two mutually antipodal 
points. 



(iii) In particular, are there manifolds with dX which is not connected? 

(iv) In particular, are there manifolds with subperiodic billiards and (or) 
with subperiodic boundary trajectories? 



3 Branching Hamiltonian flow with 
"scattering" 

Let us consider two manifolds Xi and X2 glued along the connected com- 
ponent Y of the boundary. We consider Schrodinger operators Aj on Xj 
and these operators are intertwined through the boundary conditions. We 
assume that Hamiltonian flows have simple (no branching) reflections on 
each of them, so branching comes from "reflection-refraction" . Further we 
are interested in the case when Hamiltonian flow on Xi is periodic while on 
X2 majority of trajectories are not periodic. 
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3.1 Analysis in X2 

Since we do not know any other examples of manifolds with periodic Hamil- 
tonian flows with reflections but hemisphere or half-space and A/ = 2 then 
we assume that 

(3.1) For z e T*Y ^i,r{z) = g{z) ^\ g is an antipodal map on T*Y. 

We also assume that ct>i^ and ct>2,T commute 

(3.2) cD^ ^ o 02,^ = 02,^ o cDi ^. 

Then for points of Z2,r periodicity properties of the branching Hamiltonian 
flow with reflections U^t coincide with those of ^2.t- 

Now we need to impose condition to the boundary operator except that 
{A, B} is self-adjoint. Namely, consider a point z = (x',^') G T*Y, and 
consider an auxiliary problem 

(3.3) 3j{x'.D,.e)uj = j = 1.2. 

(3.4) 5b(x',Di,e')" = 

(where u = {ui, U2)) and consider solutions which are combinations of expo- 
nents e''^^^^ where for eachj either Im Ay > or Im Aj > and {aj,xi}{x', Xj, C,') > 
0. Denote by A the set of points for which such non-trivial solutions exist. 
We assume that 

(3.5) mesr.y A = 0. 

Then we arrive to the following 

Proposition 3.1. Let in the described setup conditions (3.1), (3.2) and 
(3.5) be fulfilled. 

Let Qi, Q2 be h-pseudo- differential operators with the symbols supported 
in T*X2 and 

(3.6) |j,2,T-(n2,r n supp gi n suppq'2) = 

where 112, r is the set of point {y,ri) G ^2,t periodic with respect to ^2,t- 

Then for r((?ixe(., ., r) '^Qiy) the standard Weyl two-term asymptotics 
holds with the remainder estimate o{h^~'^). 

^^ We define 't>jr in the following way: we lift z e T*Y to z £ T^\y H Tjr with 
{aj,xi}{z) > 0, launch trajectory forward until it hits the boundary at (y', 77) G T*Xj\y 
and then project to T* Y. 
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Remark 3.2. (i) Note that 

(3.7) M,-,r(A,>) = 

where Aj ,- = Aj^ U A"^, Aj^ is the set of dead-end points of biUiard (not 
generahzed biUiard) flow and N-^ is the set of points {y.rf) E Tj^r such that 
^J^jAY'V) e ^3-j.T where Ek_-r = Lk{{x,^) E Z^,^, {a/,,Xi} = 0. 

(ii) Note that the grows properties of the branching Hamiltonian flow with 
reflections U^t coincide with those of ^2,t as long as we ensure that along 
billiards we do not approach to glancing rays on Xi: 

(3.8) {ai, xi}(x, > h'^ V(x, E 6r'^2^2,t(y, v) H ^i,. 

y{y,ri) E suppq2 ^t : ±t E [0, T] M/2,t(y, ^) e T*X2\y 

where T = T{h). Then in assumptions of section 7.4 one can recover 
similar results: namely that two-term Weyl asymptotics holds with the 
remainder estimate 0{T{h)^^h^^'^) which is usually 0(/j"'^^'^| log /)|^-^) or 
Q^^i-d+5j (depending on the growth and non-periodicity conditions to ^2,f 
Exact statement and proof are left to the reader (which is a relatively easy 
problem) . 

(iii) In (ii) we need also to assume a more sharp version of (3.5): namely if we 
consider approximate solutions to (3.3)-(3.4) with precision e we get set A,- ^ 
and we need to impose condition mesA,-;; = o{e^) or mesA^e = o{\ log£:|~-'^) 
as e — )• +0. 

(iv) In subsubsection "7.4.3" of subsection 7.4.3 we did not consider points 
z E T*Y where some of operators aj are elliptic on l^^z rather than 
hyperbolic. If we want to consider such points we need to assume (3.5) or 
more sharp version of it described in (iii). 

(v) (3.8) is guaranteed if either 32 is cylindrically symmetric on T*X2 or if 
there no internal reflection for billiards from X2 at all: 

(3.9) |{ai,xi}(x,0l>e|{a2,xi}(x,r^)|' 

V(x,0 E T*X,\y n Zi,, V(x,77) E 7*X2|y n Z2,r : 6i(x,0 = t2{x.v). 
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3.2 Analysis in Xi 

However if operators Qi and (?2 have symbols supported in Xi situation 
is very different: trajectories of ai are periodic but the typical closed 
branching billiard is the one which does not contain segments of Hamiltonian 
trajectories of 32. More precisely, due to conditions 

(3.10) ^2,r(n2,0 = 

and (3.7) we conclude that 

(3.11) ^lA^'l.r) = 

where W-^ ^ is the set of points z G 5Ii,^ such that there exists a billiard 
trajectory of U^t starting and ending in z and which is not entirely billiard 
of VUi.t. 

Then for arbitrarily large 7" > and arbitrarily small e > there exists 
a set W = Ut,s C T*X such that 

(3.12) ^ii,,(Zi,AW)<£, 

(3.13) dist(W, n;^7-UAi,,,7-) >7 

(3.14) d\st{Y,7iM)>l 

with 7 = 7(7", e) > were H'^ ^ 7- s the set of points z G Ti_r such that there 
exists a billiard trajectory of U^t of the length not exceeding T, starting and 
ending in z and which is not entirely billiard of U^i f and /\i,t,t is the set 
of points z G '^i_r such that billiard trajectory of ai hits the boundary at 
point of A2, r,r- As usual, these sets increase as T increases, their unions 
(with respect to 7~) coincide with fl'^ ^ and Ai^ respectively and Ai,^ 7 and 
^i.T.T ^^'it T ^^6 closed sets. 

Therefore there exists operator Q = Qje such that supp q C W and 

(3.15) |r((?i,e(., ., r) \{l - (?)(?2)y) - ^^'o^-' - /^'i/j'^l < Cosh^-' 

with k'. = k'jq^q^q{t). 

This estimate holds for any fixed r satisfying (3.10). Therefore we need 
to consider r((5ixe(., ., r) ^Qiy) with supp 92 C U. 

We can assume without any loss of the generality that (3.13), (3.14) are 
fulfilled for all r. Then for e < |t| < 7" 

(3.16) r{Q,,U{.,.,t)'Q2y) = V{Q,,Z{.,.,t)'Q2y) 
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where as ±t > Z(t) is Schwartz kernel of "approximate propagation with 
scattering semigroup" Z(t) constructed in the foUowing way (for a sake of 
simphcity in notations we consider t > 0: 

(i) We represent Q2 = (?2,i + ■ • • + (?2,a/ where symbols Q2,u have small 
supports; 

(ii) For each u we select = to,y < ■ • • < tM,u = t in such a way that for 
any z G supp q2,u ^t'{z) is disjoint from Y as t' = tj^^, and ^t'{z) hits Y no 
more than once as tj> < t' < tj+i^^; 

(in) We set 

Z(t)(?2,.= ( n ^e''-'^'^-^-^-'^'^^'^''^) ■ Q2,. 

0<j<M-l 

where ip G ^°° is supported in {x G Xi,dist(x, Y) > 7'} and equals 1 in 
{xGXi,dist(x, y) >2y}. 
Then 

(3.17) Z(ti + t2)Q2 = Z(t2)Z(ti)(?2 as ± ti > 0, ±t2 > 
and 

(3.18) Z{ty = Z{-t) 

and due to periodicity of U^i.t with period Tq Z{±Tq) are /)-pseudo-differential 
operators. 

Remark 3.3. In contrast to what we had before they are not necessarily 
unitary because reflection coefficients are not unitary anymore. More 
precisely, 

(i) If M^it hits T*X\y at points {x',^) with {x',^') ^ t25l2,r reflection coeffi- 
cient >iu{x',^',T) still satisfies \>iii{x' , ^' , t)\ = 1. 

(ii) However if M/it hits T*X\y at points (x',,^) with (x',,^') G i2^2,r 
then reflection-refraction matrix >c{x',^',t) = (^>cjk{x' , ^' , r)) is uni- 

tary but reflection coefficient xii{x',^',t) satisfies only |xii(x', ^', r)| < 1 
and the equality means exactly that Xi2(x', ^',t) = X2i(x', C,' ,t) = (then 
|x22(x', ^', r)| = 1 as well). 

So the principal symbol of Z(To) is the product of e'^''^^'^) where io{x,0 
is given by the usual formula for manifolds without boundary (but with 
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trajectories replaced by billiards) and also of reflection coefficients coefficients 
in the point of reflection. If the reflection coefficients do not vanish we can 
rewrite the answer in the form 

(3.19) Z{To) = e-''"^ Z{~eh-'To) = e-''* 

with /?-pseudo-differential operator L, L + L* < and e = h (but in calcula- 
tions below we are not assuming this). 

We will use this formula even if the coefficients vanish: in this case we 
simply formally allow Im £ = +oo. 

Then as x G ^o°°[e - 1, 1 - e] 

(3.20) F,^,-i,xro(t - nTo)r{Q,,Z{., ., t) Uy) = 

with 

(3.21) J{h,e,n,T) = {2nY-'^ f e'''"'«"^"^-")"+"''"^l"l) c/|Xi,,. 

We can see easily that 

(3.22) J{h, e, n, r) = o(l) as n — )■ oo provided 

(3.23) |x,({(x,0 G^r : \m£{x,^) = V Re£{x,^) = O}) = 0. 

On the other hand, note that 

(3.24) V e-^(-l"l« = ^'%'7"-'^"' as /3 > 
^ ' ^-^ l-2e-^cosa + e-2/5 

nez\o 

and formally we can extend this for /3 = as well. 

Then after summation with respect to n for /3 = we get exactly d(^[pi) 
with T defined by (6.2.63) and for /3 > we get that the expression (3.24) 
is equal to dfi'X *« {lT:)-^G{a, /3) where G{a, /3) = 2/3 {a^ + /3^y\ 

Note that 
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(3.25) Function 

(3.26) T{a,P):=T*^G{a,P) 

is harmonic function as /3 > 0, coincides with T{a) as /3 = 0, is o(l) as 

P — )■ +00 and is periodic with respect to a. 

Taking in account results of this subsection we arrive to 

Theorem 3.4. Let in the described setup conditions (3.1), (3.2), (3.5), 
(3.10) and (3.23) be fulfilled. 
Then asymptotics 

(3.27) re(.,.,r) = /s:o(r)r^+ 

(«i(r)+ / T{h-\eRei-T),-h-^€\m£)dyLr)d\ir]+o{h^~'') 

holds. 

Problem 3.5. Calculate all coefficients Kj^ for two Schrodinger operators Ai 
and A2. 

We will solve this problem in the very special case: 

Example 3.6. Assume that (at the given point z E T*Y) aj = cJ\C,\'^ and 
the boundary condition is 

(3.28) U2 = aui, 

(3.29) D1U2 = -l3DiUi 

(recall that xi > in both Xi, X2; in more standard notations of xi < in X2 
one needs to skip "— "). Then {A, B} to be self-adjoint requires ajS"^ = c^C2~^. 
We do not consider gliding points since their measure is 0. 

(i) Consider point C,' with c^l^'P < t < c||,^'p. One can prove easily that 

(3.30) >iii{e.r) = e^'^ 
with 

(3.31) v^(e',r) = arctan(|/3|-2c2-^c2(c||ef-r)^(r-c2|e'|)~^ 
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and therefore assumption (3.23) is fulfilled provided Ci 7^ C2, Cj,a,f3 are 
symmetric (i.e. CjQ = q) and the subprincipal symbol is 0. If Ci > C2 this 
case does not appear. 

(ii) Consider point ^' with cj\^'\^ < r for both J = 1,2. One can prove easily 
that 

(3.32) xn(e',r)= (cj-l)(a; + l)"' 
with 

(3.33) c.(e',r) = |/3rVq^(r-c,^|er)^(r-q^|e1)-^ 

and |xii| < 1 so assumption (3.23) is fulfilled. 

The following problems seem to be rather straightforward but not easy 
and worth publishing: 

Problem 3.7. Under stronger conditions to ^2,t and i (see (2.15)) prove 
remainder estimate 0{h^^^^^). 

Problem 3.8. Prove the same results as Hamiltonian billiards of a^ are 
assumed to be periodic only on one energy level r rather than on neighboring 
levels. 



4 Two periodic flows 

Now we consider the case of both flows M^it and ^2t generated by f{3i) and 
f (32) being periodic and satisfying (3.1) as in the following figures: 



4.1 Examples and discussion 

Let us start from an example: 

Example 4-1- Consider two Schrodinger operators in M'^ = IR+ x M.'^^^ with 
the principal symbols 

(4.1) aj{x.O = lcoj{\e + \^\'\-Ej) 

intertwined through boundary conditions. Then Tj = 2'kijjJ . 
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(a) Hemisphere, projected to equatorial plane 



(b) half-plane 



Figure 2: x-projection of billiard trajectory on two half-planes, glued 
together: original (red), reflected (blue), refracted (magenta) 



We are interested in the number N^(0) of the negative eigenvalues of 
operator As- In order to consider the eigenvalue counting function N^(A) 
for arbitrary spectral parameter A one needs to redefine 



(4.2) 



Ej 



2ujfX. 



But N~(0) is also the maximal dimension of the negative subspace of 
the operator A which does not change if we replace A by J^^AJ^^ with 
the positive self-adjoint operator J and the latter equals to N^(0) for the 
latter operator. Picking up J = J{x) equal coj in Xj we arrive to the case of 
ojj = 1 and 7j = 27r. Therefore in this example one can make periods equal 
T"i = 7"2. The dependence oi uj will come through Ej redefined by (4.2). 

There are however many deficiencies of the described approach; most 
important is that it does not work with the following example (unless some 
unnatural conditions to coj, Ej are imposed): 

Example J^.2. Consider two operators on hemisphere §^ with the principal 
symbols 



(4.3) 



^i(^.0 



UJ; 



Ej) 
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where |^| is calculated in the standard Riemannian metrics. 

The main obstacle here is that multiplying operator by J^z from both 
sides and calculating f{A) do not play well together. 

4.2 Reduction to the boundary 

First of all notice that the contribution to the remainder of the domain 
{z, dist(z, T*Y) < Eq} does not exceed Coelh^'^ and therefore one needs to 
consider the contribution of {z, dist(z, T*Y) > sq} with arbitrarily small 
but fixed Eq. 

In this zone one can consider U{x,y, t) as the solution of the Cauchy 
problem with respect to xi (or yi) with data at {xi = 0} (or {yi = 0}). 
Rewriting therefore U = R^U^Ry where R is an operator resolving this 
Cauchy problem is given by an oscillatory integral we can rewrite 

(4.4) V{U{., ., t) %) = \-'{Q{x'. hD'^. hD„ h)V) 

with /)-pseudo-differential operator R and V = Sx,m-iSy,m-i^ the Cauchy 
data for U (assuming that m is an order of A). 

One can see easily that as Q is an operator with the symbol supported 
in T*Xj, the principal symbol of Q at (x', ^, r) is defined as an averaging of 
Q along ^jt{zj) with t E [0, tj{zj)] where zj = r'^{x',^') n Tj^ and tj{zj) is 
the time of the next hit of the boundary. 

Then instead of continuous family of Fourier integral operators e''' ^'^ 
we can consider a discrete family J-"" defined by the following way: consider 
Vj = dm-iUj = (i/j~, Vj^) with Vj^ corresponding to incoming and outgoing 
solutions. 

We have two types of trajectories: those in Z^v which hit T*X\y at 
points elliptic for {bs-j — t) where _/' = 1, 2 and those which hit T*X\y at 
points hyperbolic for [a^ — r) for both k = 1,2. The analysis of the former 
is of no different from what we have seen in the previous subsection: we just 
consider U^yt and construct the (real-valued) symbol i as long as we assume 
that (x', ^) ^ A where A is defined in the paragraph preceding (3.5). 

The analysis of the latter is more interesting. Therefore as before we 
have a matrix of reflection-refraction {>Cjk)j,k^i,2 which is a symbol defined 
on W X [t — e' , T + e') where W is a a zone described above. This matrix is 



Chapter 4. Two periodic Sows 20 

unitary in the norm 

(4.5) \\V\\ = {J2{\v;f + \Vrf)\{aj,x,}U^of ■ 

J 

We also have diagonal unitary matrices M^ = diag(e'^^'", e'^^") with iji, 
calculated along z/-th leg of the closed trajectory of ^jt with t = | Tq and 
i/ = 1,2. 

Consider closed branching billiard originated from point z E T*Y; first 
there are two trajectories (in Xi and X2) both hitting T*X\y at g{z) and 
then there are two trajectories (in Xi and X2) returning to z. Let us recall 
that g{z) is an antipodal point. 

Therefore branching billiard is characterized by an unitary (in (4.5)-norm) 
matrix 

(4.6) M{z, r) = x(z, r) M'(z, r) x(^(z), r) M'{g{z), r). 

4.3 Analysis of the evolution 

One can prove easily that modulo Ojih^^'^) 

(4.7) F,^,-i,XT{tr{Q{x', hDl, hD„ h)V) ^ 

A//ij=o-i,...j„)e{i,2}" 

where x ^ "^0°° i^ supported in [|, 1] and symbols Mjk and Q and Tj are 
calculated at {x',^',t) where Tj = Tj{t) since all trajectories are periodic 
on energy levels close to r. 

One can rewrite terms with equal N as 

(4.8) (2nhf-' I x{\)e'''-''^^n^-^^^.) I I ^ 

i=0i,-J/v)e{i,2}" 

Mj,j,Mj^, . . . Mj,_j,e-''-'^^''n+-^^^)Qj^j^dx'd^'dX 
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which in turn equals 

(4.9) {27r hf-' I x{\) II tr{S{\,T,x',a^Q)dx'd^'d\ 
where 

(4.10) 5{XT .-.-.0=(^,-(Ar-.--)r.^^^ e'(----^)-=m J 

and M = {rriji^). 

Consider eigenvalues of this matrix; as 7 — )• oo they tend to eigenvalues 
of the same matrix with /i = 0. If 

(4.11) mes{(x,^') : p e Spec{5{0,T,x ,^'))} = o{l) 

Vr G [r - e/?T"\r + e/?7"^] as T -)■ cx) 

then this term is o{h^~'^) and therefore we estimated expression (4.7) by 
o{Th^~^) which in the end of the day returns remainder estimate o{h^~'^) 
with the Tauberian main part. 

Note that instead of S we can consider matrix 

(4.12) S{r,x',C)=[ _,,,4* 

\^ m2i e g22^ 

with 7* = |(7i-72) which is different from S(0, ., ., .) by factor e2'''"^^(''"i+'^2)_ 
Consider correction to the main part; we are interested at r = 0; so 
we integrate (4.7) with x = 1 from r = — ootor = resulting in (after 
multiplication by h~^ and modulo o(h^^^)) 



(4.13) {2nhr'xT{T, + ---+T,) J2 E 

^ -^ A//ije{i,2}" 

Mjj,Mj^s ■ ■ ■ Mj,_j,{Tj, + ... TjJ-^Qj,j,dx'dC 

where we plug r = 0. 
This latter equals to 

(4.14) Q(r) = {27rhy-' II tr T(/.(x', ^' , h-\)) Q{x' , ^.T)dx'd^' 
where L is defined by e''- = S and T is defined by (6.2.63). 
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Theorem 4.3. Consider two Schrodinger operators Ai and A2 with all 
periodic trajectories on levels close to t, satisfying (3.1). Further, let (3.5) 
and (4.11) be fulfilled. Then asymptotics 

(4.15) tr r(e(., ., r)) = Xo(r)/)-^ + (xi(r) + n{r))h'-' + o{h'-') 
holds. 

Example 4-4- Consider example 3.6 in the new conditions to Haniiltonian 
flows. We can use resuhs of that example immediately to treat billiards 
with complete internal reflections (thus non-branching). 

To treat branching billiards we need to calculate eigenvalues of the matrix 
(xy/f). It follows from calculations of example 3.6 that X22 = — Xn; thus we 
arrive to 

(4.16) A = e^"^, (y9 = arccosxii = arccos((cj — l)(a; + l)^"*"; 
then obviously (4.11) is fulfllled and asymptotics (4.15) holds. 

Remark 4-5. On the contrary assume that (4.11) is not fulfilled. Then there 
will be eigenvalues of high multiplicity or clusters of eigenvalues located in 
o(/))-vicinities of solutions r of equation det(S — 1) = 0. One can rewrite 
this equation as 

(4.17) cos(-v9+ -r/?"^(7i+ T2)) = cosa-cos(-^+ -r/)"^(7i - T2)) 

e'('t>+^)cosa e''('^+^) sin a 
-e'^-^-^^^sina e''^"^^^) cos a 
unitary matrix. 



as M = { ^,.(^_^) .(^_^) ^ ) which is the general form of the 



The following problem seems to be a difficult one and worth of publica- 
tion: 

Problem 4-6. Prove the same results as fiows ^j_t are assumed to be periodic 
only on one energy level r rather than on neighboring levels. 

On the contrary, the following problem seems to be relatively easy: 

Problem 4- 7. (i) Prove the same results as there are more than two manifolds 
Xj j = 1, ... , m' with the all billiards periodic. 

(ii) Extend these results to the case when there are also manifolds Xj 
j = m' + 1, ... , m with almost all billiards non-periodic. 
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